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5 METHOD OF VISUALIZING SOUND FIELDS OF INDIVIDUAL SOUND 

SOURCES USING ACOUSTIC HOLOGRAPHY 

Priority Claim 

[1] This application claims priority from Korean Patent Application No. 10-2003- 

10 0002354 filed 14 January 2003, which is herein incorporated by reference. 

Field of the Invention 

[2] The present invention relates generally to a method of visualizing the sound 

fields of sound sources using acoustic holography, and more particularly, to a method which 
discriminates the sound fields of individual sound sources having the same frequency 
15 component and visualizes the sound fields thereof. 

Description of the Related Art 

[3] Generally, acoustic holography is a method of predicting acoustic quantities, 

such as sound pressure, sound velocity and sound intensity, at arbitrary positions by 
measuring sound pressures on a two-dimensional plane called a hologram plane. The 
20 positions and radiation characteristics of sound sources can be detected using such 
acoustic holography. The acoustic holography has been widely researched and used in 
various fields, such as the automobile industry, munitions industry, and educational 
laboratories for researching acoustics. 

[4] Near-field acoustic holography allowing actual sound sources to be detected 

25 was proposed by Williams and Maynard, and is disclosed in a thesis, "Holographic imaging 
without the wavelength resolution limit" of Phys, Rev. Let., 45, 554-557, 1980 by E.G. 
Williams and J.D. Maynard, and research on the application and improvement thereof have 
continuously progressed. 
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5 [5] Several conventional methods are disclosed below for methods of measuring 

sound pressures on a hologram plane. First, a method of simultaneously measuring sound 
pressures at all measurement points is disclosed in U.S. Pat. No. 4,415,966 entitled 
"Nonwavelength-limited holographic sound field reconstruction" issued in 1983 to J. D. 
Maynard and E. G. Williams. Second, a method of sequentially measuring sound 

10 pressures is disclosed in a thesis, "STSF-a unique technique for scan-based near-field 
acoustic holography without restrictions on coherence" of B&K Technical Review No. 1 , 
1989, by J. Hald. Third, a method of continuously measuring sound pressures is disclosed 
in a thesis, "Moving frame technique for planar acoustic holography" of J. Acoust. Soc. AM., 
103, 1734-1741, 1998 by H. S Kwon and Y. H. Kim. Fourth, a system and method for 

15 visualizing acoustic characteristics by measuring hologram of a moving sound source is 
disclosed in Korean Pat. No. 217872 issued in 1998 to Y. H. Kim and H. S. Kwon. 

[6] Further, a system and method for visualizing band-limited noise and a system 

and method for monitoring the state of transportation machinery using the same are 
disclosed in Korean Pat. Appl. No. 366206. 

20 [7] An acoustic holography method measures sound pressures using planar 

acoustic holography, cylindrical acoustic holography, spherical acoustic holography and the 
like according to the shape of a hologram plane, and is disclosed in a thesis, "Nearfield 
acoustic holography: I. Theory of generalized holography and development of NAH" of J. 
Acoust. Soc. Am., 78, 1395-1413, 1985 by J. D. Maynard, E. G. Williams and Y. Lee. 

25 [8] Further, an acoustic holography method, which can be applied to a hologram 

plane with arbitrarily shaped surfaces on the basis of a boundary element method, has 
been developed and is disclosed in a thesis, "Digital holographic reconstruction of sources 
with arbitrarily shaped surfaces" of J. Acoust. Soc. Am., 85, 588-598, 1 989 by W. A. 
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5 Veronesi and J. D. Maynard. 

[9] The above research on the acoustic holography mainly focuses on obtaining all 

sound fields. If there are multiple sound sources having the same frequency component, 
only the sum of individual sound fields other than the individual sound fields of individual 
sound sources can be obtained based upon the research. 

10 [10] Meanwhile, if multiple sound sources exist, conventional methods of obtaining 
sound fields created by individual sound sources under the condition where a sound field is 
stationary and sound sources are incoherent, are described below. 

[11] Hald attempted to obtain individual sound fields using singular value 
decomposition, which is disclosed in a thesis, "STSF-a unique technique for scan-based 

15 near-field acoustic holography without restrictions on coherence" of B&K Technical Review 
No. 1, 1989 by J. Hald. Further, Hallman and Bolton attempted to obtain individual sound 
fields on the basis of a partial coherence method, which is disclosed in a thesis, "Multi- 
reference near-field acoustic holography" of Proc. Of Inter-Noise 92, 1 165-1 170, 1992 by D. 
Hallman and J. S. Bolton. Further, Kwon and Bolton proposed a method of obtaining 

20 individual sound fields using a residual method, which is disclosed in a thesis, "Partial field 
decomposition in nearfield acoustical holography by the use of singular value 
decomposition and partial coherence procedures" of Proc. Of Noise-Con 98, 649-654, 1989 
by H. S. Kwon and J. S. Bolton. Ruhala and Burroughs proposed a method of correcting 
the partial coherence method, which is disclosed in a thesis, "Separation of leading edge, 

25 trailing edge, and sidewall sound sources from rolling tires" of Proc. Of Noise-Con 98, 109- 
114, 1998 by R. J. Ruhala and C. B. Burroughs. 

[12] However, these methods are problematic in that sensors must be arranged 
close to the individual sound sources to detect in advance the positions of the sound 
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5 sources. Further, these methods are problematic in that, since the positions of the sound 
sources are detected in advance before the acoustic holography is performed, these 
methods are contradictory to the acoustic holography aiming to locate the positions of 
sound sources. 

[13] A method of obtaining sound fields created by individual sound sources in a 
10 frequency domain under the condition where a sound field is stationary and the sound 
sources are incoherent will be described below. 

[14] Fig. 1 is a block diagram of a typical single-input/single-output system, and Fig. 2 
is a block diagram of a typical two-input/single-output system. 

[15] In Fig. 1 , an input signal X(f) is distorted due to a noise component N(f) and then 
15 output as an output signal 9(f) . At this time, the noise component N(f) is incoherent to the 
input signal X(f). If the output signal 9(f) and the input signal X(f) are simultaneously 
measured, the auto-spectrum of an output signal Y(f), not distorted, can be obtained by 
Equation [1], 

^ xx U ) 

20 where Sxx(f) is the auto-spectrum of the input signal X(f), and Sxv(f) is the cross spectrum of 
the input signal X(f) and the measured output signal 9(f) . 

[16] If a signal linear or coherent to the input signal X(f), such as R(f) of the following 
Equation [2], instead of the input signal X(f), is measured and then applied to Equation [1], 
Equation [1] can be expressed by Equation [3]. 

25 R(f) = C(f)X(f) [2] 

^ RR U ) 
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5 [17] In Fig. 2, if two input signals Xi(f) and X 2 (f) are independently operated, that is, 
incoherent, the auto-spectrum of an output signal Y(f) is expressed by the sum of the 
contributions of the individual input signals or sound sources as shown in Equation [4], 

(/) [4] 

where Syiyi (f) is the contribution of the first input signal, and Sv2Y2(f) is the contribution of 
10 the second input signal. 

[18] Since the X^f) and X 2 (f) are incoherent, the output obtained by the X 2 (f) is only 
noise from the point of view of the Xi(f). Therefore, if the output signal Y(f) and the first input 
signal Xi(f) or a signal Ri(f) coherent to the Xi(f) are measured, the contribution of the first 
input signal Xi(f) can be obtained by Equation [5]. 

m ~s rn~~s 77V 1 J 

[19] As described above, if multiple independent or incoherent sound sources exist in 
an acoustic system, input signals of individual sound sources or signals coherent to the 
sound sources must be obtained so as to obtain the contributions of individual sound 
sources to the output spectrum. 

20 [20] The input signals of individual sound sources or signals coherent thereto can be 
obtained from vibrations or sound pressure signals generated close to the individual sound 
sources. As described above, the method of obtaining sound fields created by the 
individual sound sources using acoustic holography is disadvantageous in that the sound 
fields can be obtained only when the vibrations or sound pressure signals of the individual 

25 sound sources are measured while sensors such as reference microphones are arranged 
close to the individual sound sources. 

[21] With all the above descriptions taken into consideration, the conventional 

5 
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5 acoustic holography method can obtain only the sum of sound fields of the individual sound 
sources other than the individual sound fields thereof if multiple sound sources having the 
same frequency component exist. 

[22] Moreover, the conventional acoustic holography method is problematic in that, 
since it obtains individual sound fields using signals measured through sensors arranged 
10 close to the sound sources, the positions of the sound sources must be detected in 
advance before holography is performed. 

Summary Of The Invention 

[23] Accordingly, the present invention has been made keeping in mind the above 
problems occurring in the prior art, and an object of the present invention is to provide a 

15 method of visualizing the sound fields of individual sound sources using acoustic 

holography, which evaluates the value of a maximum sound pressure on a sound source 
plane as a signal coherent to a specific sound source using sound pressures measured on 
a hologram plane, and obtains the sound field of the sound source using the coherent 
signal without measuring the positions of the sound sources before acoustic holography is 

20 performed. 

[24] In order to accomplish the above object, the present invention provides a method 
of visualizing sound fields of individual sound sources using acoustic holography, 
comprising the steps of a) calculating sound pressures on a sound source plane using 
sound pressures measured on a hologram plane; b) extracting a maximum sound pressure 
25 on the sound source plane, and evaluating a value of the sound pressure at the position as 
a signal coherent to a specific sound source; c) obtaining a sound field of the extracted 
sound source using the coherent signal; d) eliminating the sound field of the extracted 
sound source from sound fields of all sound sources, and determining whether any 
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5 remaining sound field exists; and e) returning to the step b) if any remaining sound field 
exists at step d). 

[25] Preferably, the sound field of the extracted sound source at step c) is obtained by 
calculating contribution of each of the sound sources using the coherent signal in a spectral 
matrix consisting of auto-spectra at individual points on the sound source plane and cross 
10 spectra between different points. 

Brief Description Of The Drawings 

[26] The above and other objects, features and other advantages of the present 
invention will be more clearly understood from the following detailed description taken in 
conjunction with the accompanying drawings, in which: 

15 [27] Fig. 1 is a block diagram of a typical single-input/single-output system; 

[28] Fig. 2 is a block diagram of a typical two-input/single-output system; 

[29] Figs. 3a and 3b are a conceptual view and a flowchart of a signal processing 
method to obtain sound fields of individual sound sources according to the present 
invention, respectively; 

20 [30] Fig. 4a is a perspective view of an acoustic holography system according to an 
embodiment of the present invention; 

[31] Figs. 4b and 4c are upper and side views of a speaker experiment device for 
verifying the acoustic holography system of Fig. 4a, respectively; 

[32] Fig. 5 is a graph showing the spectrum of a sound pressure measured by a 
25 single reference microphone if only a single speaker is driven in Fig. 4a; 

[33] Fig. 6 is a distribution chart of the sound fields of sound pressures on a sound 
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5 source plane of Fig. 5, obtained by a conventional method; and 

[34] Figs. 7a to 76 are distribution charts of the sound fields of sound pressures of 
individual speakers of Fig. 5, obtained by the present invention. 

Description Of The Preferred Embodiments 

[35] Hereinafter, a method of visualizing sound fields of individual sound sources 
10 using acoustic holography according to embodiments of the present invention will be 
described in detail with reference to the attached drawings. 

[36] Figs. 3a and 3b are a conceptual view and a flowchart of a signal processing 
method to obtain the sound fields of individual sound sources according to the present 
invention, respectively. 

15 [37] Referring to Fig. 3b, the sound pressures on a hologram plane are measured 
(step S300). 

[38] Next, a total sound field of a two-dimensional sound source plane on which 
sound sources are arranged is obtained using an acoustic holography method (step S310). 

[39] Thereafter, the position where a sound pressure has a maximum value on the 
20 sound source plane is determined as the position of a first sound source, and the value of 
the sound pressure of the first sound source is selected as a signal coherent to the first 
sound source (steps S320 and S330). 

[40] Further, a sound field created by the first sound source is obtained using the 
coherent signal selected at step S320 (step S340). 

25 [41] Further, it is determined whether any remaining sound field exists after the sound 
field created by the first sound source is eliminated from the total sound field (step S350). 
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5 [42] Thereafter, if any remaining sound field exists after the first sound field is 
eliminated from the total sound field, the remaining sound field steps S320 to S350 are 
repeatedly applied to (step S360). 

[43] The above process is repeatedly carried out until all of the individual sound fields 
created by the individual sound sources are obtained. 

10 [44] The above-described process is carried out so that, after the individual sound 
sources and signals coherent thereto are selected from a spectral matrix including of auto- 
spectra at individual points on the sound source plane and cross spectra between different 
points, the sound fields created by the individual sound sources are discriminated from each 
other by using the coherent signals. 

15 [45] Hereinafter, the above process is described in detail with reference to equations. 

[46] If sound pressures at N positions on the sound source plane, measured at step 
S300, are Psi(f)> Ps2(f),-, PsN(f), respectively, a spectral matrix S S s(f) consisting of auto- 
spectra at individual points on the sound source plane and cross spectra between different 
points is expressed by the following Equation [6], 
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where E is an expected value, and * is a conjugate complex number or a conjugate 
transpose operation. 

[47] At step S31 0, the performance of acoustic holography at a stationary random 
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5 sound field is equal to the obtainment of the spectral matrix as shown in Equation [6] for 
obtaining sound pressures or sound velocities on the sound source plane or other planes 
by measuring sound pressures on the hologram plane. 

[48] Sss(f) of Equation [6] can be expressed by the sum of the contributions of M 
sound sources as shown in the following Equation [7], similar to Equation [4] obtained in the 
10 case where M independent or incoherent sound sources exist, 

°ss (f) [7] 

where Ssisi(f) is the contribution of an i-th sound source. 

[49] A rank of S S s(f) of Equation [6] is determined by the number of independent 
sound sources, and the value thereof is M. Further, since S S s(f) is Hermitian and positive 
15 semi-definite, it can be decomposed as the following Equation [8], 

s ss (f) = r;(f)r s {f) [8] 

where r 5 (/) is a M x N matrix, and is not uniquely determined. The Sss(f) expressed by 
Equation [7] can be decomposed using any decomposition method, such as a singular 
value decomposition or LU (Lower triangular matrix and Upper triangular matrix) 
20 decomposition method. 

[50] In Equation [8], if an m-th column of the r s (/) is expressed by r s m (/) , the auto- 
spectrum of a sound pressure at m-th position on the sound source plane can be 
expressed by , m (f)r Sm (f) . Further, if the position where the auto-spectrum has a 
maximum value is ow, and an m ma x-th column of the r s (f) is r S max (/) , the contribution of 
25 the first sound source to the Sss(f) can be obtained by the following Equation [9]. 

£ (y) _ s, m ax (/)r 5 max (fWs (f) ^ 

* 6, max 
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5 [51] In the spectral matrix Sss(f) consisting of auto-spectrums at individual points on 
the sound source plane and cross spectrums between different points, a spectral matrix 
Sss-i(f) consisting of remaining contributions after the contribution of the first sound source 
expressed by Equation [9] is eliminated can be obtained by Equation [10]. 

10 [52] As described above, a procedure for obtaining a spectral matrix consisting of 
remaining contributions after the contribution of one sound source is eliminated from the 
spectral matrix Sss(f) consisting auto-spectra at individual points on the sound source plane 
and cross spectra between different points, and a procedure for eliminating the contribution 
of another sound source from the spectral matrix and obtaining a spectral matrix consisting 

15 of the remaining contributions, are repeatedly performed. Therefore, the contributions or 
sound fields of the individual sound sources can be obtained. 

[53] Fig. 4a is a perspective view of an acoustic holography system according to 
an embodiment of the present invention, and Figs. 4b and 4c are upper and side views of 
a speaker experiment device for verifying the acoustic holography system of Fig. 4a, 
20 respectively. 

[54] Referring to Figs. 4a to 4c, six speakers 1 0 to 60 are arranged within an acrylic 
box 100, and operated independently. Further, in order to measure spectrums for 
frequency bands of sounds radiated from the individual speakers 10 to 60 and provide 
reference sound pressures, six reference microphones (not shown) are arranged to be 
25 spaced apart from the center of the acrylic box 1 00 by approximately 1 m. Further, in order 
to measure the pressures of sounds radiated from the individual speakers 10 to 60, a 
microphone array 70 is provided. 

[55] Further, in order to process attribute data of sound sources, measured through 
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5 the reference microphones and the microphone array 70 and visualize a hologram, a signal 
processing means, such as a computer 80, is provided. 

[56] In the acoustic holography system having the above construction, a hologram 
plane is positioned to be spaced apart from the center of the acrylic box 1 00 by 42.7cm, and 
a sound source plane is positioned to be spaced apart from the center of the acrylic box by 
10 32.1cm. 

[57] After the acoustic holography system is constructed as described above, sound 
pressures on the hologram plane are measured in a vertical direction using the microphone 
array 70 in which 24 microphones are arranged at intervals of 5 cm in the vertical direction. 
Further, the sound pressures on the hologram are measured at 60 positions in a 
15 circumferential direction while rotating the microphone array 70 by stages at an angle of 6 
degrees around the acrylic box 100. 

[58] Fig. 5 is a graph showing the spectrum of a sound pressure measured by a 
single reference microphone if only a single speaker is driven in Fig. 4a. 

[59] Fig. 6 is a distribution chart of the sound fields of sound pressures on the sound 
20 source plane of Fig. 5, obtained by a conventional acoustic holography method. 

[60] Referring to Fig. 6, a frequency is 400Hz, and four speakers are primary sound 
sources in this frequency band. However, as shown in Fig. 6, only sound fields created by 
three speakers are vaguely shown and it is difficult to observe a sound field created by the 
4th speaker. Further, it can be seen that even the sound fields of the three speakers, 
25 observed vaguely, are not clearly discriminated from each other due to a spatial overlap 
phenomenon. 

[61] Figs. 7a to 7d are distribution charts of the sound fields of sound pressures of the 
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5 individual speakers of Fig. 5, obtained by the present invention, and respectively show the 
distribution charts of the sound fields created by the first to fourth speakers. 

[62] Referring to Figs. 7a to 7d, even a sound field of a speaker, which could not be 
observed in Fig. 6, can be observed, and the sound fields of other three speakers can be 
more clearly observed. 

10 [63] As described above, the present invention provides a method of visualizing 
sound fields of individual sound sources using acoustic holography, which discriminates 
individual sound sources from each other and visualizes the sound fields thereof using 
sound pressures of the sound sources, obtained using acoustic holography, even though 
the positions of the sound sources are not known. Therefore, the present invention is 

15 advantageous in that there is no need to provide sensors close to the individual sound 
sources, the resolution of acoustic holography can be improved, and the priorities of the 
individual sound sources can be determined. 

[64] Further, the present invention is advantageous in that, even if the positions of the 
sound sources have been detected in advance, that is, even in the case where it is difficult to 
20 arrange sensors close to the sound sources as in the case of the measurement of automobile 
tire and engine noises and aero-acoustics, the present invention can discriminate individual 
sound sources from each other by applying a software approach. 

[65] Although the preferred embodiments of the present invention have been disclosed 
for illustrative purposes, those skilled in the art will appreciate that various modifications, 
25 additions and substitutions are possible, without departing from the scope and spirit of the 
invention as disclosed in the accompanying claims. 
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